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Abstract 



We perform next-to-leading order global analyses of deep inelastic and related data 
for different fixed values of as(Mg). We present sets of parton distributions for six values 
of as in the range 0.105 to 0.130. We display the (x,Q 2 ) domains with the largest 
parton uncertainty and we discuss how forthcoming data may be able to improve the 
determination of the parton densities. 



Analysis of the scaling violations of deep inelastic scattering data provides one of the accu- 
rate ways to determine the QCD coupling as- The scaling violations observed in recent high 
precision muon and neutrino deep inelastic data yield values[] of as = 0.113 ± 0.005 |[J and 
0.111 ± 0.006 respectively. Moreover a global parton analysis which includes these data, 
with other related data, gives as = 0.113 ± 0.005 ||. However, there are other independent 
determinations which lie outside this range; for example as determined from LEP event shapes 
or from r decays. A recent review of all the methods is given by Webber ||J] who concludes 
that the world average is 0.117±0.005. 

as is now also being determined from jet rates observed in the experiments at HERA |J 
and Fermilab |J. These methods have the advantage of determining as over a wide range of 
Q 2 and it is believed that they will eventually yield values with equal precision to the other 
determinations of as- However, they require the use of parton distributions which have their 
own particular value of as and so the question of consistency arises. Does the "output" as 
depend on the "input" value of «s? In order that the sensitivity to the input partons may 
be studied, we repeat the global analysis of refs. || [7|] for various fixed values of as (Mf ) in 
the interval which covers these other independent determinations, namely 0.105 to 0.130. This 
study allows us to highlight the deep inelastic data that particularly constrain a s . Moreover 
it provides a quantitative estimate of the uncertainty associated with the parton distributions 
fi(x, Q 2 ) in different regions of x and Q 2 . 

An earlier analysis ||, which studied the uncertainty in the determination of as, did provide 
4 parton sets which cover a limited range of as- Here we extend the range and use the improved 
set of deep inelastic and related data. CTEQ || have recently presented an additional parton 



set with a high as, and Vogt (10| has provided 5 sets of GRV |Tl[ partons with as in the interval 
(0.104, 0.122). The latter is not a global analysis and so cannot accommodate the variation 
of partons, particularly at larger x, which attempt to compensate for the shift of as from its 
optimum value. 

We base our global next-to-leading analysis on the MRS (A) parametric forms, updated in 
ref. J7j to include recent HERA data [|12|, That is we consider variations about the optimum 
fit MRS(A')Q Fig. 1 shows the x 2 values for various subsets of deep inelastic data obtained 
in six new global fits^ with different values of as(M§). Cross section and asymmetry data 
for Drell-Yan and W hadroproduction are included in the analysis, but their contributions to 
X 2 are not shown. These data remain well- described in the fits with different as values by 
slight adjustments of the partonic flavour structure of the proton; they do not pin down as- 
Also the data for prompt photon production are accommodated as as varies by a change in 

throughout we work in terms of the QCD coupling at the Z pole, Q 2 — M§, evaluated in the MS scheme 
with 5 flavours, which we shall denote simply by as- 

2 We choose not to base our analysis on the MRS(G) set of partons [(7), since this would involve an extra 
parameter which is only relevant to data at very small x. It therefore would introduce a degree of ambiguity in 
a domain where the precision of the data is going to rapidly improve. We return to this point later. 

3 The six sets have as = 0.105,0.110, . . .0.130 and are denoted by MRS. 105 etc. The parton sets can be 
obtained by electronic mail from W.J.Stirling@durham.ac.uk. 
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scale. The x 2 profile for the WA70 data is shown in Fig. 1. Only the smallest values of as 
are disfavoured; the prediction does not fall off quite steeply enough, as the photon transverse 
momentum increases, to agree with the observed distribution. 



The HERA 1993 data |L2| |13| are included in the analysis. In addition we show the x 
values for the preliminary 1994 ZEUS data that were obtained with the electron-proton 
collision point shifted (so that the detectors can gather deep inelastic events at smaller x). Due 
to the logarithmic scale that has been used for x 2 it is easy to be misled by Fig. 1 about the 
relative importance of various data sets in the determination of as- The x 2 profiles at the top 
of the plot have a more significant impact than those which lie lower down. 

Considerable insight into the effect of varying as (and the related ambiguities) can be 
obtained from Fig. 2. This shows the available data for F 2 P = F 2 P at three particular x 
values: x = 0.0008 in the HERA range, x = 0.05 which is relevant for W production at 
Fermilab and x = 0.35 representative of the large x BCDMS precision data which provide the 
tightest constraints on as- The curves are obtained from the three parton sets which have 
as = 0.105, 0.115 and 0.125. Recall that the optimum overall description occurs for as = 0.113 
and so the continuous curve gives a better global fit than the ones either side. As expected, the 
scaling violation is greatest for the partons with the largest value of as- Also, as may perhaps 
be anticipated, the curves cross in the region of the data, which lie in different intervals of Q 2 
for the different values of x. Away from the (x, Q 2 ) domain of the data the predictions show a 
considerable spread. For example for x = 0.0008 and Q 2 ~ 10 3 GeV 2 we see quite a variation 
in the prediction for F 2 P . The ambiguity in the small x domain is actually greater than that 
shown, since the quark sea and the gluon have been constrained to have the same small x 
behaviour, that is 

xS ~ A s x~ Xs , xg ~ A g x~ Xs (1) 

with = X g . Unfortunately HERA is unable to measure F 2 in this region of x and Q 2 ; the reach 
of the collider is kinematically limited to the domain x/Q 2 > 10~ 5 GeV -2 . Nevertheless as the 
precision of the HERA data improves it will be possible to determine As and X g independently 
(see 0). The sensitivity of the predictions to the interplay between the form of the gluon and 
the value of as demonstrates the importance of a global analysis which includes the crucial 
large x constraints on as- The x 2 profiles shown in Fig. 1 that are obtained from the HERA 
data overconstrain as since they are based on fits which set As = A s . In our global "as" 
analysis this has a negligible effect on the partons, except at small x, where for Q 2 values away 
from the HERA data there will be more variation than the spread that our curves imply. 



The lower plot in Fig. 2 shows a typical set of the high precision BCDMS data ||14[| . In the 
large x domain these data place tight constraints on as, free from the ambiguity associated 
with the gluon. 

Fig. 3 gives another view of the constraints on the partons in various regions of x and Q 2 . 
It is a contour plot of the ratio R of F 2 as predicted by two sets of partons with as(M|) fixed 
either side of the optimum value. In (x, Q 2 ) regions of precise data, acceptable fits demand 
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that the ratio R be equal to 1. This is strikingly borne out. We see that the R = 1 contour 
lies precisely in the centre of the band of the fixed-target deep inelastic scattering data. Again 
R ~ 1 in the region of the more accurate HERA data, that is x < and Q 2 ~ 15 GeV 2 . 
When the precision of the HERA data improves and the accuracy extends over a larger domain 
of x and Q 2 we would expect the R ~ 1 contour to also track the HERA band, but probably at 
the expense of allowing and X g to be free independent parameters. Of course Fig. 3 is just 
an overview of the description of one observable (F 2 ep = F^). The global fit has many other 
constraints to satisfy. Nevertheless _F 2 ep is measured over far wider regions of x and Q 2 than 
the other observables and so Fig. 3 gives a useful indication of features of the global fit. 

Additional experimental data in regions where the contours in Fig. 3 are closely spaced would 
clearly have significant impact on the analysis. For example sufficiently precise information in 
the region of high x and low Q 2 (the lower left corner of Fig. 3) would help pin down as even 
further. This is the domain of the SLAC experiments ||2(J . However we must take care to avoid 
regions where there are appreciable target-mass/higher- twist effects. For the SLAC data these 
effects are smallest in the region x ~ 0.3, provided that Q 2 > 5 GeV 2 @, 0, ||. The Table 
below shows the x 2 values for the subset of the SLAC data |20) that lie in the "safe" region 
(0.18 < x < 0.45) obtained from our sets of partons with different as'- 



a s (M 2 ) 


0.105 


0.110 


0.115 


0.120 


0.125 


0.130 


X 2 (25 pts) 


37 


24 


25 


46 


93 


179 



The SLAC data clearly support the optimum value of as determined by the other deep inelastic 
data, that is as = 0.113. The inclusion of the SLAC data in the global fits would simply 
mean that the curves at x = 0.35, for example, would cross-over at a value just below the 
Q 2 = 50 GeV 2 intersection shown in Fig. 2. 

Returning to Fig. 3 we see that jet production at large Et at Fermilab samples partons 
in an x, Q 2 domain far removed from the regions where deep inelastic measurements exist. 
For example, jets produced centrally with a transverse energy Et = 300 GeV sample x = 
2Et/\/~s ~ 0.3 and Q 2 ~ E T - Fig. 4 compares the observed single-jet inclusive distribution of 



CDF [[23[] with the predictions from partons with three different values of as- For simplicity 
we have evaluated the leading-order expression at a common renormalization and factorization 
scale /i = Et/2. Of course a precision comparison between data and theory will require a full 
next-to-leading order analysis p4j| . However, our aim here is to compare the spread of the 
predictions with the uncertainty of the data. A leading-order calculation is entirely sufficient 
for this purpose. A change of scale simply boosts the predictions up or down relative to the 
data, but leaves the shapes in Fig. 4 essentially unchanged. To consider the implications for 
partons it is necessary to discuss the description in two distinct regions of Et- For Et ^ 200 
GeV the jet cross section is dominated by the gg and qg initiated subprocesses. Here the cross 
section ratios reflect the different shapes of the gluon distributions in the region 0.05 < x < 0.2, 
with the predictions spread out even more by the differences in the associated values of a 2 s . 
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For Ex ^ 200 GeV, on the other hand, the cross section becomes increasingly dominated by 
quark-initiated subprocesses. Here the differences between the curves reflect the spread in the 
predictions of F 2 at large x (x > 0.25) and large Q 2 (Q 2 ~ E T ), but now suppressed (rather than 
enhanced) by the differences in a 2 s {Q 2 ). This illustrative exercise demonstrates the value of a 
precise measurement of the jet distribution. If the experimental uncertainties can be reduced 
then these data will impose valuable constraints on the gluons at small x (x ~ 0.1) and on the 
quarks at large x (x ~ 0.35), as well as on the value of as- 

In summary we note that deep inelastic scattering data determine as to be 0.113 ± 0.005. 
This value is found in an analysis of the BCDMS and SLAC data by Milsztajn and Virchaux 
[[L| and in the global analyses |3], [7| which include, besides the BCDMS measurements, other 



deep inelastic and related data. Moreover the SLAC deep inelastic data [20], which are not 
included in the global fit, also favour this value of as, see the Table above and the sample data 
in Fig. 2. The deep inelastic determination of as relies mainly on the scaling violations of the 
data in the large x domain (x > 0.2), a region free from the gluon and its ambiguities. It is 
easy to verify that the low x HERA deep inelastic data for do not determine as unless 
restrictive assumptions about the small x behaviour of the gluon and sea quark distributions 
are made. Indeed we find that the HERA data give little constraint on as even with the 
assumption that As = X g in (p. The values of as determined from other (non deep- inelastic) 
processes cover a wider interval: 0.110 < «g < 0.125 [|J. Some methods rely on input partons 
and so there is a need for parton sets with values of as which cover this interval. We have 
therefore performed a series of global analyses of the deep inelastic data to obtain realistic sets 
of partons corresponding to a sequence of values of o;s(M|). Since as is not optimal these are 
compromise fits, which yield curves that intersect in the x, Q 2 regions where precise data exist; 
see, for example, Fig. 2. The body of the deep inelastic data occurs in the region Q 2 ~ 20 GeV 2 . 
Fig. 5 shows the gluon and up quark distributions from 3 parton sets with different as for Q 2 
values above and below this value. The systematics displayed in these plots may be anticipated 
from Fig. 2. For low Q 2 , below the body of the data, we see from Fig. 5(a) that the lower as 
partons "swing more about the x ~ 0.05 — 0.1 pivot points in favour of lower x" ; and vice- versa 
for the high Q 2 partons of Fig. 5(b). It is interesting to note that W and Z boson production 
at the Fermilab pp collider sample u and d partons with x ~ 0.05 and Q 2 ~ 10 4 GeV 2 , that is in 
the region of the pivot point. The predicted production cross sections are therefore unusually 
stable to the change of the set of partons used in the calculation^. Jet production, on the other 
hand, samples partons over a range of x and Q 2 , as well as being directly dependent on as(Q 2 ). 
It is therefore important to have to hand sets of realistic partons with different as(M|). 
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Figure Captions 



Fig. 1 



Fig. 2 



Fig. 3 



Fig. 4 



The x 2 values versus the value of as used in the global analysis. The contributions to \ 2 
are shown for the BCDMS 0, CCFR 0], NMC WA70 0, HI (1993) and 
ZEUS (1993) data sets. The x 2 values for the preliminary ZEUS (SVX) 1994 JTj 



data are also shown, but these data are not included in the fit. The x 2 values are also 
shown for the MRS(A') set of partons [0], which has the optimum value of as, namely 
a s = 0.113. 

The scaling violations of _F 2 ep = F^ v at three different values of x. The data are from 
refs. [0, [T(| [T2|, |T3|, [TS], fH| . The curves correspond to the global parton fits with as = 
0.105, 0.115 and 0.125. At x — 0.35 some data points not normally included in our global 
fits are also shown: first, the lower Q 2 BCDMS measurements which are made only at 
their lower beam energy and, second, SLAC data |2(J in the region which is insensitive 
to target- mass/higher-twist corrections. 

Contours of fixed R = F 2 ep (as = 0.125)/F 2 ep (as = 0.105) in the x, Q 2 plane, where 
^2 P ( a s) is the structure function calculated from partons obtained in a global analysis 
with as fixed at the given value. The bands indicate the regions where measurements of 
F 2 exist. The HERA data are much more precise towards the low Q 2 end of the band. 

The pp-initiated jet Et distribution at y/s = 1.8 TeV normalized to the prediction from 
partons with as = 0.115 (i.e. MRS. 115). The data are the CDF measurements of 
d 2 <j I 'dETdrj averaged over the rapidity interval 0.1 < \rj\ < 0.7 f23fl . The curves are 
obtained from a leading-order calculation evaluated at rj = 0.4. The data are preliminary 
and only the statistical errors are shown. The systematic errors are approximately 25% 
and are correlated between different Et points. We thank the CDF collaboration for 
permission to show these data. 



Fig. 5 The xg and xu parton distributions at (a) Q 2 
parton sets with as = 0.105, 0.115 and 0.125. 



5 GeV 2 and (b) Q 2 = 10 4 GeV 2 of the 
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